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ABSTRACT 

A large sample of young stellar groups is analysed aiming to investigate their cluster¬ 
ing properties and dynamical evolution. A comparison of the Q statistical parameter, 
measured for the clusters, with the fractal dimension estimated for the projected clouds 
shows that 52% of the sample has substructures and tends to follow the theoretically 
expected relation between clusters and clouds, according to calculations for artificial 
distribution of points. The fractal statistics was also compared to structural parame¬ 
ters revealing that clusters having radial density profile show a trend of parameter s 
increasing with mean surface stellar density. The core radius of the sample, as a func¬ 
tion of age, follows a distribution similar to that observed in stellar groups of Milky 
Way and other galaxies. They also have dynamical age, indicated by their crossing 
time, that is similar to unbound associations. The statistical analysis allowed us to sep¬ 
arate the sample into two groups showing different clustering characteristics. However, 
they have the same dynamical evolution, since the whole sample has been revealed as 
expanding objects, for which the substructures seem to have not been erased. These 
results are in agreement with simulations that adopt low surface densities and models 
under supervirial conditions. 
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1 INTRODUCTION 

It is widely accepted that the physical conditions for the 
formation and evolution of stellar clusters can be inferred 
on basis of embedded clusters structure, which depends on 
the distribution of dense gas of their original cloud. Numer- 
ous clusters e mbedded in clouds are found in the Galaxy. 
lLada fc Ladal J2003h suggested that most of them shall lose 
their dynamical equilibrium, been dissolved in a few tens of 
Myr. Based on the surface density profile of the distribution 
of members, they proposed two types of cluster structure: 
centrally concentrate d and hi erarch ical. 

lElmegreen fc Falgaronei dl996l f demonstrated the frac¬ 
tal origin of the relations of size and mass distribution in 
the interstellar gas, determined from cloud s surv eys in the 
literature. More recently, iGirichidis et al.i (120121 ') used hy¬ 
drodynamic simulations of formation of filaments and pro¬ 
tostars, in different initial configurations, aiming to evaluate 
the statistical properties of embedded clusters in the early 
star-forming phases. They found that more substructures 
are formed by flat density profiles and compressive modes, 
instead of to produce centrally concentrated profiles. 

Classical open clusters are characterized by a centrally 
concentrated surface distribution, whose radial profile is 
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smooth and can be reproduced by a simple power law func¬ 
tion. On the other hand, hierarchical type clusters show 
multiple peaks, or substruct ures in their surface density dis¬ 
tribution, like NGC 1333 ( Ladaj^lves&Lada Il996h . for 
instance. More recently, iGutermuth et ahT i 2008 ) discussed 
the structure of this cluster based on Spitzer observations 
of 137 members, whose surface density map shows an elon¬ 
gated distribution, not centrally concentrated, suggesting a 
low velocity dispersion and that NGC 1333 is not in virial 
equilibrium. 


In spite of the large number of different clusters found in 
the Galaxy, their relative occurrence remains to be verified, 
requiring the investigation of clustering properties in a as 
large as possible sample of studied objects. 

For instance, iFeigelson et al.l (1201 il l studied the struc¬ 
ture of a huge sample of young stars detected by ray-X ob¬ 
servations in the Chandra Carina Complex Project (CCCP). 
They found three different types of clustering properties in 
their sample that includes previously known clusters and 
some additional smaller groups identified by CCCP. The 
classical, centrally concentrated clusters, present an uni- 
modal stellar distribution. Objects showing several substruc¬ 
tures, distributed in a circular boundary, are classified as 
multimodal clumpy structure, which is suggested to be an 
unequilibrated stage of cluster formation. The third type of 
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cluster contains sparse compact groups, without any con¬ 
centration that is possibly due to triggered star formation. 

The effects of ionizing fronts on triggered star forma- 
tion h ave be en ev aluated using numerical simulations by 
IWalch et al.l (|2013l ). for instance. 

Considering the rapid changes on the spatial distri¬ 
bution of forming stars shown by these simulations, it 
is expected a strong evolution of the clustering struc- 
ture, t h at can b e eva luat ed by the statistical parameter Q 
dCartwright fe Whitworthl 120041) . Indeed, a rapid change of 
Q is noted in_the simulat ions o f i onization- induced star for¬ 
mation bv lDale. Ercolano fe Bonnelll (j2012i . l2013i 'l. However, 
smooth evolution is also found in the simulations with- and 
without fee dback of ionizing sou r ce, depending on th e initial 
conditions (IParker fe Daldl2013l : IParker et all 2014T ) . 

Based on fractal statistics, in a previous work we dis¬ 
cussed the formati on and evolution con diti ons of four young 
clusters (iFernandes, Gregorio-Hetem fc Heteml 120121 . here¬ 
after FGH12). The preliminary results revealed that one ob¬ 
ject of the sample has radial stellar density profile that is not 
similar to the fractal characteristics of the nearby cloud, in 
a contrary way of the other three studied clusters. 

Aiming to validate, refine and enlarge our previous in¬ 
vestigation of fractal statistics and its possible relation to 
the initial conditions of cluster formation, in the present 
work the analysis is extended to a sample of 21 open 
star clusters, selected according to their youth and in¬ 
termediate distances. The clust ers w ere characterized by 
ISantos-Silva fe Gregorio-Heteml d2012l . hereafter SG12) by 
using near-infrared data to derive their structural param¬ 
eters based on stellar density maps , similar to the study 
developed bv IBonatto fe Bical (|2009|) . for instance. 

The dependence of fractal statistics on age and 
crossing time was also analysed, aiming to investi¬ 
gate the dynamical evolution that is expected for the 
star distribution in young clusters, which has been 
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among others). The results wer e also compared 
to th e observatio na l stu di es by Schmeia fe Klessenl 
1 20061). ICaballerd (120 0^1 . Sanchez fe Alfarol (120091) . 
iPortegies Zwart. McMillan fe Gielesl (l2010l . hereafter 
PZMG10) and lWright et ali ( 2014f) '. for instance. 

The paper is organized as follows. Section[2]is dedicated 
to present the characteristics of the selected clusters. Sec¬ 
tion [3] summarizes the adopted methodology, describing the 
calculations of fractal dimension of clouds, clustering statis¬ 
tics of surface stellar density, and error estimation by using 
the bootstrap method. Section [3] presents the results from 
the statistical analysis, and in Sect. [5] the dynamical evolu¬ 
tion of the sample is discussed. Finally, the main conclusions 
are summarized in Sect.[S] 


2 SELECTED STAR CLUSTERS 

All the selected young clusters have similar angular sizes 
(R < 20 arcmin), and intermediate distances, which is d ^ 
2 kpc for most of the objects. Since the observed E(B-V) is 


Table 1. Sample studied in the present work and informations 
from literature. 


Cluster 

N'j' 

Age 

d 

n 

R 

r c 



(Myr) 

(pc) 

(pc- 2 ) 

(pc) 

(pc) 

Collinder 205 

174 

5 

1800 

4.8±1.9 

3.4=1=1.2 

0.30=1=0.03 

Hogg 10 

88 

3 

2200 

7.6±1.7 

1.9=1=0.4 

0.46=1=0.28 

Hogg 22 

97 

3 

1700 

9.0±3.8 

1.9=1=0.5 

1.9=1=1.1 

Lynga 14 

68 

3 

950 

16.7=1=7.0 

1.1±0.4 

0.22=1=0.04 

Markarian 38 

57 

6 

1500 

6.8±2.3 

1.6=1=0.5 

0.13=1=0.03 

NGC 2302 

70 

3 

1700 

4.Oil.4 

2.3=1=0.5 

0.54=1=0.10 

NGC 2362 

124 

4 

1480 

3.3=1=1.0 

1.9=1=0.4 

0.26=1=0.05 

NGC 2367 

60 

3 

2200 

2.3=1=0.6 

2.9=1=0.7 

0.34=1=0.11 

NGC 2645 

104 

5 

1800 

7.5=1=2.1 

2.1±0.5 

0.22=1=0.04 

NGC 2659 

215 

5 

2000 

5.2=1=0.9 

3.6=1=0.5 

1.93±0.37 

NGC 3572 

149 

3 

1900 

7.7=1=1.9 

2.5=1=0.5 

0.16=1=0.03 

NGC 3590 

79 

3 

1680 

11.7=1=4.5 

1.5=1=0.4 

0.29=1=0.06 

NGC 5606 

98 

3 

2200 

4.8=1=0.9 

2.6=1=0.4 

0.51=1=0.13 

NGC 6178 

106 

5 

1430 

12.2=1=2.7 

1.7=1=0.4 

0.30=1=0.12 

NGC 6604 

90 

6 

1600 

4.4=1=1.3 

2.6=1=0.7 

0.50=1=0.11 

NGC 6613 

133 

5 

1550 

4.9=1=1.2 

2.9=1=0.7 

0.12=1=0.01 

Ruprecht 79 

174 

5 

2700 

3.6=1=0.5 

3.9=1=0.5 

2.20=1=0.80 

Stock 13 

65 

4 

2000 

3.8=1=0.9 

2.3=1=0.3 

0.11=1=0.05 

Stock 16 

139 

7 

2000 

5.8=1=1.6 

2.8=1=0.7 

1.20=1=0.60 

Trumpler 18 

164 

5 

2850 

2.5=1=0.4 

4.6=1=0.8 

0.38=1=0.10 

Trumpler 28 

73 

2 

1050 

12.3=1=4.5 

1.4=1=0.5 

0.84=1=0.16 

Berkeley 86 

85 

5 

1585 

3.4 

2.8 


NGC 2244 

295 

3 

1660 

1.9 

7.0 


NGC 2264 

292 

3 

760 

5.0 

4.3 


NGC 6530 

62 

3 

1300 

2.9 

2.6 



Note: Data for 21 clusters are given by SG12, and for the last 
four clusters are from FGH12 and references therein. 


low, varying from 0.1 to 0.7, we assume that none of them 
can be considered an embedded cluster. 

Tabic Q] lists the 25 young star clusters analysed in the 
present work, giving the information available in the litera¬ 
ture. For comparison purposes, the four clusters previously 
studied by FGH12 have been included in the sample. The 
other 21 clusters have structural and fundamental parame¬ 
ters determined by SG12. 


2.1 Observed radial density profile 

SG12 used surface stellar density to derive the structural 
parameters of the clusters. Members hip probabi lities, evalu¬ 
ated according to proper movement (IDias et al.ll2002l) , were 
used to identify the very-likely members. SG12 also adopted 
color-magnitude filters that enable to separate field-stars 
from candidate members, which were also included in the 
total number of objects ( Nt ) used to determine the stellar 
density maps. 

Table [l] reproduces the structural parameters obtained 
by SG12 by fitting the observe d rad ial density profile, based 
on the empirical model from iKingj d 1962|) : a(r) = a bg + 
00/(1 + (r/r c ) 2 ), where cr(r) is the stellar surface density 
(stars/arcmin 2 ), a bg is the average density measured in a 
reference region (background), no and r c are respectively 
the density and the radius of the cluster core. The cluster 
radius ( R ) is given by the point where <r(r) = <j bg . The 
adopted errors are the uncertainties provided by the fitting 
method. In Sect. [4] these parameters are compared to the 
statistical fractal parameters obtained in the present work. 
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2.2 Mass and age 

By comparing the clu s ter members with pre-MS 
( Siess. Dufour fe Forestinil 120001 ') and main sequence 
( Girardi et al.l 120021 ') models in the unreddened near-IR 
colour-magnitude diagrams, SG12 estimated fundamental 
parameters like age and mass. The mean ages are ~ 5 Myr, 
but some of the clusters have a bimodal age distribution, 
showing a second peak around 12 Myr (Hogg 22, Lynga 14, 
NGC 2302, NGC 2645, NGC 5606, and NGC 5606). 

The observed mass distribution of each cluster was fit¬ 
t ed w ith the mass function oc from iKroupal 

(l200lll . The x values determined by the fitting method in¬ 
dicate that 33% of the sample have low power-law ind ex, 
comparable to the slope \ = 0.3 ± 0.5 assumed bv iKroupal 
d200lf) for M < 0.5Mq. The other 67% of the clusters have 
deeper slopes x ~ 0-8 1° 1-5 that are, within the errors, 
roughly similar t o the Sa l peter ’s IMF (y = 1.35), which 
is the index that iKroupal (1200 ill suggested for M > 1M 0 . 
This result is also in agreement with the mass distribution 
with x ~ 1.0 that has been fou nd for the initial cluster mas s 
function of diverse clusters (e.g. lElmegreenll 20 ld : IOevll 201 lh . 

The mean mass of the cluster members is > 1 M@, in¬ 
dicating an incompleteness of low-mass stars, probably due 
to the 2MASS detection limit that constrains the presence 
of faint sources in the sample. 

Considering the total mass, it is interesting to note that 
our sample corresponds to an inte rmediate range bet ween 
the embedded clusters studied bv lLada fe Ladal l)2003h and 
the massive “l eak” clusters h aving R > 10 pc in the sam¬ 
ple studied bv lPfalzneil (l201lll . All of them follow the same 
mass-radius relation given by M ~ 118i? 1,3 . SG12 (see their 
Fig. 4c) co mpare the sampl e with the distribution of clusters 
analysed bv IPfalznerl (| 201 lh who proposed that these objects 
represent the ending of a time sequence, which mass-radius 
dependence is M ~ 359.R 1 ' 1 . However, this power law fit¬ 
ting is significantly different of the distribution mentioned 
above, found by SG1 2 for our samp l e, as well as for the em¬ 
bedded clusters from lLada fe Ladal d2003l l and also for clus¬ 
ters sh owing R < 10 pc in the sample analysed bv IPfalznerl 
ll 201 lh . 

The results f rom SG12 fa vour the int erpretation given 
by lAdams et all (l2006l l and lAd am 1 dm3) suggesting that 
differences on mass-radius relations are more likely due to 
different formation conditions instead of a time sequence 
endin g in the massive “leak” clusters proposed bv IPfalzneit 

l 201 lh . 

A similar conclusion was achieved by FG12 from the es¬ 
timated volumetric density, which is correlated with cluster 
radius by p = 28 R~ 17 , where p is given in M 0 /pc 3 * * * * , which 
stands for ou r samp le as well as for embedded clusters from 

lLada fe Ladal d2003ll . be ing quite similar to other s tudies , 

like p oc R~ 192 found bv ICamargo, Bonatto fe Bical l2Qlolj . 
for instance. 


3 METHODOLOGY 

FGH12 performed fractal analysis of the clusters Berke¬ 

ley 86 , NGC 2244, NGC 2264 and NGC 6530. By comparing 

with the statistic derived from artificial data simulating stars 

distributions and parental clouds, the clustering parameters 


of three of the analysed clusters are compatible with the 
fractal dimension of the observed clouds. NGC 6530 is the 
only cluster showing a central concentration of stars that 
is different of the substructures found for the nearby cloud. 
In the present work our previous analysis is extended to a 
larger sample, following the same methodology adopted in 
FGH12, which is summarized in this section. 


3.1 Fractal dimension of the nearby clouds 

Visual extinction (Ay) maps were used to perform statistical 
analysis of the clouds. Aiming to investigate the interstel¬ 
lar matter distribution in the direction of the clusters, the 
Av maps we re extracted fr om the Dark Clouds Catalogue 
presented bv lDobashi et all |2005h . 

For each cluster we searched for 2° x 2° images contain¬ 
ing the Av levels obtained from 2MASS data, with spacial 
resolution of ~ T per pixel. A threshold of S/N = 10 was 
chosen to define the minimum level of Ay contour, avoiding 
low density regions in the calculations. Therefore, Av maps 
are useful to indicate the presence of dense clouds nearby 
the cluster. 

Figure |T] shows some examples of Av maps and cor¬ 
responding spatial distribution of cluster members, which 
were used in t he fractal analysis discus sed in Sect. 13.21 

Following iHetem fe Lepinel dl993l l. the fractal dimen¬ 
sion D 2 of contour levels is measured by using the perimeter- 
area relation of the clouds P = (A 1 ^ 2 / K)° 2 , where P is the 
perimeter of a given Av contour level and A is the area 
inside it, and K is a parameter connected to the geome¬ 
try (shape) of the objects under study. The uncertainties on 
the estimates of D 2 and K of projected images of clouds 
were o btaine d by applying standard least squares methods 
dPress et al . IM). 

Since this method uses bi-dimensional projected maps, 
the comparison of D 2 with the values of fractal dimension 
from literature (D 3 ), determined from three-dimensional 
distributions, requires a 3 D x 2D t ransfo rming function 
(D 3 = _F(Z? 2 )). Following iFalconeil (|l990t l. we adopt the 
relation D 3 = D 2 + 1 , which is a good ap pr oxima tion for 
D 3 < 2.5 dCartwright, Whitworth fe Nutteil l2006h . How¬ 
ever, the relationship between fractal dimension of a cloud 
projection and three dimensional calculations may depend 
on a more comple x inference of D 3 = T(D 2 ) th at is not 
clearly established dSanchez, Alfaro fe Perezll2005l l. 


3.2 Clustering statistics 


The parameter Q, which carries information about the 
cluster fractal structure, is a dimensionless quantity given 
by Q = rn/s, where m and s are statistical parameters 
that d epen d _on the geometric distribution of data points 
dGartwright fc Whitworthll2004 ‘'l. 

The parameter to is the normalized mean edge length, 
related to the surface density of cluster members projected 
position, defined by: 


JV-l 


(AvIV )!/ 2 

2=1 


where N is the total number of considered points, rrii is the 
length of edge i in the minimum spanning tree (described 
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Figure 1. Left: Visual extinction (Ay) map (galactic coordinates) of nearby clouds of four objects of the sample. Contours represent 
Ay levels with S/N > 10. The location of the star clusters is indicated by the central rectangle. Middle: A zoom of the Ay map showing 
the spatial distribution of the cluster members. Right: Minimal spanning tree and smallest circle (dashed line) that contains all members 
of the cluster. Upper panels show clusters with Q > 0.8 (NGC 6530 and NGC 2367), while small scale substructures (NGC 2244 and 
NGC 2264) are seen in lower panels. 
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below) and An is the area of the smallest circle that contains 
all points projected on the plane of the cluster. Each smallest 
circl e was d e termi ned by adopting the algorithm proposed 
by iMegiddol (1 19831 1 . 

The minimal spanning tree is defined as the unique 
network of straight lines that can connect a set of points 
without closed loops, such that the sum of all the lengths 
of these lines (or edges) is the minimal possible. Figure [Q 
gives examples of smallest circle and respective spanning 
tree, which wer e cons tructed using the method described by 
iGower fe Rossi dl969li in order to obtain m for each cluster. 

The value of s represents the mean separation of the 
points, and is given by 

N -1 N 

s= N{N-1 )R n 5Z l^-^l 

i=l j=l + i 

where n is the vector position of point i, and Rn is the 
radius of the smallest circle that contains all points. 

Studies of the hierarchical structure in young clusters 
have used the Q paramet er to distingui sh fragmented from 
smooth distributions (e.s. lElmegreenfeoift l. Comparing the 
cluster physical structure with the geometry of its possible 
parental cloud may tell us about the original gas distribution 
of the star forming region and how the particular cloud may 
have evolved during the cluster formation. 

Using _ artific ial distri bu tions of points, 

iLomax, Whitworth fe C artwright) (|20 1lf) adop ted the 
technique proposed by ICartwright fc Whitworth! d2004) to 
perform a statistical analysis that gives inferences on the 
fractal dimension measured on grey-scale images generated 
by models of stars clusters and clouds. They used data sets 
varying from 64 to 65536 points, which give a parameters 
space of Q as a function of radial profiles (index a) and 
fractal dimension ( D 3 ) that are simila r to, but larger than 
the c alculati ons performed b y ICartwright fe Whitworthl 
(|2004l f and [Sanchez fc Alfarol (|2009l 4 . for instance. In 
Sect. 14.21 the results are compared with these theoretical 
models. 

3.3 Standard deviation estimated with 
bootstrapping 

The bootstrap method was adopted in order to estimate the 
standard deviation of the statistical param eters. This nu¬ 
merical method, introduced bv lEfroril (jl979h . is independent 
of model or calculations; is not based on asymptotic results; 
and is simple enough to be automatized as an algorithm. 
The method creates a set of simulated data by considering 
each parameter of the model as varying in a chosen con- 
hdence level according to a given distribution. Then, this 
artificial set is treated as a real sample (as those obtained 
by observations) and the standard deviation is calculated 
from these data. 

In the present work, the b ootstrap pi ng alg orithm was 
implemented as proposed bv IPress et alj (119921 1. The algo¬ 
rithm generates synthetic distributions of cluster members 
based on the original positions of the stars. A number of 
stars is moved from their original positions, being replaced 
by random points within a box, around the star, with side 
equals to the cluster radius. Then, the parameters rn, s and 
Q are calculated for each synthetic distribution. By adopting 
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Figure 2. Distribution of Q as a function of number of stars 
(top axis) and membership probability (bottom axis) calculated 
for members and candidate members of NGC 2244. The limiting 
values adopted in the present work, N = 256 and P = 50%, are 
shown by a dashed line. 


a fraction f = 1/e ~ 37% of replaced stars in a symmetric 
normal distribution, it was possible to derive la deviation 
for Q, m and s. 

4 ANALYSIS OF FRACTAL STATISTICS 
4.1 Membership dependence 

In spite of the whole set of stars (candidates and members) 
had been used by SG12 to determine the cluster structural 
parameters, in the present work we preferred to constrain 
the statistical analysis to confirmed members only. Based 
on the probability membership (P%), the stars were sep¬ 
arated in two sub-samples: the very likely members of the 
cluster have P ^ 50% (denoted by P m ), while the candidate 
members (called P c ) have P < 50% . 

Aiming to evaluate the effect on the fractal statistics 
by considering different numbers of stars and to explore the 
validity in choosing P m to represent the very likely members, 
Fig. [5] shows the variation of Q as a function of N for the 
cluster NGC 2244. The adopted result Q = 0.74 obtained 
for N = 256 members is plotted with error bars. It can 
be seen that this result stands as a representative value for 
almost the whole data set, with a smooth decreasing for low 
numbers of stars. A significant change in Q occurs only for 
N = 14 members with P > 80%, but the main conclusion 
about the clustering structure of NGC 2244 remains the 
same. 

Figure [5] also displays the variation of Q as a function 
of the membership probability. This distribution shows that 
the choice of Pm is good enough to avoid large errors due to 
a small data set (low N), also excluding possible field-stars 
contamination (low P%). 

For each set of stars (P m and P c sub-samples), we per¬ 
formed the calculation of the parameters m and s, which 
are displayed in Fig. [3] It can be seen that the distribution 
of Pm members in this plot tend to follow the expected cor¬ 
relation between m and s, while the candidate members do 
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S 

Figure 3. Trend of m increasing with s shown by the very likely 
members of the sample (SG12 and FGH12). Candidate members 
(P c ) do not follow this trend that is illustrated by Q = 0.8 (full 
line) and 10% deviation (dashed lines). Data for other clusters 
from literature are also plotted for comparison (see Sect. HS- 
med symbols represent Q < 0 . 8 . 
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Figure 4. Parameter Q, estimated from the spatial distribution 
of cluster members, as a function of fractal dimension (-D 3 ) of 
projected clouds. Typical error bars are shown for NGC 6530. 
For comparison, clustering statistics obtained for artificial data 
(shaded areas) and other theoretical results (full lines) are also 
shown. 


not follow the same distribution, probably due to field-stars 
contamination. 

For comparison with literature, other objects are also 
included in Fig- El like p Oph, Tau rus, Chamaeleon, IC 348, 
Serpens (iSchmeia fc Klessen 20061). a nd IC 2391 which was 
analysed by Cartwright fc Whitworthl (I2004T) along with the 
first four clusters mentioned above. In spite of the agreement 
on the Q values between both works, their results on m and 
s differ in almost 20%, excepting Taurus for which the m 
values are 50% discrepant._ 

Since the cluster a Ori (ICaballerol l2008i ) has available 
parameters that can be used in several of our comparative 
plots, a different symbol is used distinguishing it of other 
clusters from literature. 


4.2 Cluster structure compared with projected 
clouds 

Table[2]gives the statistical parameters, along with the num¬ 
ber of studied stars in each cluster, and an estimation of 
galactic coordinates (l,b). 

Aiming to derive the size of the area occupied by the P m 
members, we determined the mean (l,b) that give the center 
of the stellar spatial distribution. The respective dispersion 
is used to estimate the observed sizes, which show a circular 
symmetry (Al ~ A 6 ) suggesting that our objects are not 
elongated. 

I 11 Fig. [4] the parameter Q, obtained from the Pm 
spatial distribution, is compared with the fractal dimen¬ 
sion of the nearby clouds (discussed in Sect. 13.11) . The lo¬ 
cus of statistical parameters calculated for artificial data 
points ((Lomax. Whitworth fc Cartwrightl l2011) is shown as 
a shaded area. Above Q = 0.8 (dotted line) the shaded 
area shows the dependence on radial profiles (index a), 
while the variation of fractal dimension ( D 3 ) is found be¬ 
low this line. Since the index a was not estimated for our 
objects, the upper region of Fig. [4] is displayed only for 


comp arison with other theoretical results (I Sanchez fe Alfarol 
120091 )(full lines). They are also comparable to simulations 
for cluster types F2.0 to F3.0 and 3D0 to 3D2.9 from 
ICartwright fe Whitworthl (120041 ). which respectively corre¬ 
spond to fractal and radial distributions. 

It can be noted that 43% of the sample is found around 
the bottom shaded area of Fig. [4] meaning that substruc¬ 
tures observed for these clusters (Q < 0.8) coincide with the 
fractal characteristics of the projected clouds (D 3 < 2.5). 
It is important to stress that this is not necessarily a clue 
for physical relation with the nearby cloud. The meaning 
of these clusters coinciding with the theoretically expected 
Q x D 3 parameters space is a possible similarity of the sub¬ 
structures found for the clusters and the clouds that are 
projected in their direction. 

The opposite occurs for the clusters that have Q > 0.8 
(empty circles), including NGC 6530 (indicated by error 
bars), which do not show the same clustering structure ob¬ 
served in the projected clouds. The values of Q obtained for 
these clusters suggest they are more centrally concentrated, 
while the nearby clouds have small scale substructures. 

In fact, the original geometric structure of the gas dis¬ 
tribution is expected no longer be the same of the remain¬ 
ing cloud after the cluster formation. It also depends on the 
initial conditions of the star-forming environment. Further¬ 
more, the investigation of clusters evolution requires an esti¬ 
mate of their dynamical age, in order to verify if their initial 
structure should have been erased or not (see discussion in 
Sect. [SJ. 


4.3 Surface density 

The structural parameters calculated by SG12 were com¬ 
pared to Q, but no special trend was found for most of them. 

We also searched for correlations among structural pa¬ 
rameters compared to the clustering parameters m and s in 
separate. With the exception of NGC 3572 and a Orion, a 
trend was found for clusters having Q > 0.8, which show s 
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Figure 5. Mean surface density compared to s, which varies more 
with radial clustering. Symbols are the same as Fig. [3] A trend 
of s increasing with n is observed only for clusters with Q > 0.8. 


Age [Myr] 

Figure 6. Core radius as a function of age of our clusters (SG12, 
filled symbols for Q <0.8) and a Ori, compared to Associations 
and Young Massive Clusters of Milky Way (MW), and clusters of 
other galaxies of the Local Group (LG) and outside the LG. 


increasing with average surface density n (number of stars 
per pc 2 ). The same does not occur for the other clusters that 
have substructures, as seen in the n x s plot of Fig. [5] Also, 
no trend was found in the comparisons with m. These are 
expected results, since s is more sensitive to radial cluster- 
ing, while m varies more under the presence of substructures 
dCartwright fc Whitworth]|2004l ). 

In order to investigate if clustering is related to the clus¬ 
ter concentration parameter, which is given by tidal radius 
over core radius, the results were compared to t he rel ation 
Q vs. log(R t /r c ) obtained by I Sanchez fe Alfard ll2009l) . 

According to the definition, r c is the radius at which the 
surface brightness drops by a factor of two from the central 
value (PZMG10). 

The tidal radius of th e sample was estimated b y usin g 
the relation suggested by ISaurin. Bica fe Bonattol J2012I) : 
Rt = (M/M ga i) 1/3 d G c, with M gat = Vq C <Igc/G, where 
<1gc is the distance of the cluster to th e Galactic center , 
Vgc = 254 km/s and Rgc = 8.4 kpc dReid et al.ll2009l) . 
The values of Rt are in the range of 7.6 to 13 pc, which give 
concentration parameter log(Rt/r c ) that varies from 0.7 to 
2, i ndependently of Q . 

ISanchez fe Alfaro! d2009l ) analysed a sample of Galactic 
open clusters. For the subsample with Q > 0.8 they found a 
correlation that is not confirmed by our clusters. Contrary 
to their results, the distribution of all objects in the sam¬ 
ple is not distinguished, having a behaviour similar to the 
fractal-like subsample ana lysed by them. The same is noted 
for a Ori (lCaballerdl2008l ) that has Rt =21 pc, r c =l pc and 
Q=0.88. These results possibly indicate that our sample is 
not tidally limited, meaning that tidal radius is not well de¬ 
fined for them. Yet, Rt estimates depend on many assump¬ 
tions having large error bars. For this reason, we preferred 
not using tidal radius in the present analysis. 

5 DYNAMICAL EVOLUTION 

No special trend was found when comparing statistical pa¬ 
rameters with mass and age of the clusters from SG12. This 
unexpected lack of correlation is probably due to the simi¬ 
larities observed in their fundamental parameters. 

Therefore, the analysis of age dependence on cluster 


structure is better illustrated by comparing the sample with 
objects spread in a larger range of fundamental parameters 
values, like the young massive clusters (YMC) and stellar 
associations of Milky Way and other galaxies, using data 
available in Tables 2 to 4 of PZMG10, for instance. Also in 
this section, the Q evolution and the crossing time as a func¬ 
tion of age are compared with observations and theoretical 
results from literature. 

5.1 Core radius vs Age 

It is noted in Fig.[6]that the sample has age and core radius 
consistent with the Milky Way YMCs, but a large dispersion 
is observed. The inclusion of our objects, as well as a Ori, 
in this distribution makes less clear the expected trend of r c 
increasing with age. However, our results confirm the predic¬ 
tion by PZMG10 that young clusters with large r c probably 
were missing in their sample, due to a selection bias. 

It is also interesting to note that our sample tends to fill 
the upper region in Fig. [S] (large r c values) that is the ex¬ 
pected locus of unbound associations, for which there are 
more available data of size than r c (see Fig. 8 (left) of 
PZMG10). 

The hypothesis that our clusters probably are expand¬ 
ing objects can be confirmed based on their dynamical evo¬ 
lution. In the following subsection, the size and mass esti¬ 
mated by SG12, were used to investigate how many crossing 
times the clusters have undergone, which gives an indication 
on how dynamically evolved they are. 

5.2 Crossing time 

By using the total mass (M) and cluster radius ( R ) 
from SG12, the crossing time of the clusters were esti¬ 
mated adopting the expr e ssio n T cr = 10 (R 3 /GM) 1/2 from 
iGieles fe Porteeies ZwartJ (1201 ll ) . They suggested a bound¬ 
ary distinguishing stellar groups under different dynamical 
conditions that is expressed by the “dynamical age”, also 
called parameter n, which is given by the ratio of age and 
crossing time. Unbound associations (expanding objects) 
have n < 1, while bound star clusters have n > 1. 
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Table 2. Results from fractal analysis 


cluster 

1 ° 

Al 

b° 

A b 

N 

S 

m 

s 

d 2 

Te r (Myr) 

Collinder 205 

269.21 

0.12 

-1.85 

0.11 

34 

0.84T0.12 

0.65T0.21 

0.77T0.30 

1.13L0.05 

47.4 

Hogg 10 

290.8 

0.05 

0.08 

0.05 

18 

0.80T0.12 

0.79T0.40 

0.98L0.53 

1.24L0.09 

26.5 

Hogg 22 

338.56 

0.07 

-1.14 

0.06 

23 

0.81±0.11 

0.76T0.38 

0.94T0.49 

1.38T0.03 

23.1 

Lynga 14 

340.92 

0.07 

-1.09 

0.07 

15 

0.75T0.18 

0.74T0.32 

0.99T0.53 

1.40T0.03 

15.3 

Markarian 38 

11.98 

0.06 

-0.93 

0.06 

15 

0.75T0.19 

0.67T0.22 

0.89T0.35 

1.31T0.06 

25.8 

NGC 2302 

219.3 

0.07 

-3.12 

0.07 

30 

0.76T0.15 

0.67T0.25 

0.88T0.42 

1.17L0.05 

40.2 

NGC 2362 

238.17 

0.07 

-5.55 

0.06 

22 

0.81±0.11 

0.74T0.30 

0.91L0.37 

1.34L0.07 

21.4 

NGC 2367 

235.59 

0.07 

-3.83 

0.07 

19 

0.81T0.14 

0.69T0.25 

0.84L0.36 

1.48L0.04 

59.2 

NGC 2645 

264.8 

0.07 

-2.9 

0.08 

28 

0.80T0.12 

0.72T0.29 

0.90L0.43 

1.39L0.04 

29.4 

NGC 2659 

264.18 

0.11 

-1.65 

0.11 

55 

0.83T0.12 

0.73T0.27 

0.88L0.40 

1.39L0.04 

45.8 

NGC 3572 

290.71 

0.07 

0.2 

0.07 

29 

0.83T0.13 

0.64T0.17 

0.77L0.26 

1.24L0.10 

30.6 

NGC 3590 

291.21 

0.05 

-0.17 

0.05 

13 

0.84T0.17 

0.83T0.44 

0.99L0.45 

1.24L0.09 

20.7 

NGC 5606 

314.85 

0.06 

0.99 

0.07 

32 

0.82T0.14 

0.75T0.32 

0.92L0.47 

1.36L0.04 

38.6 

NGC 6178 

338.41 

0.08 

1.21 

0.07 

21 

0.76T0.15 

0.72T0.29 

0.95L0.48 

1.34L0.03 

22.3 

NGC 6604 

18.25 

0.09 

1.7 

0.1 

33 

0.74T0.19 

0.65T0.19 

0.89T0.40 

1.42L0.03 

37.5 

NGC 6613 

14.18 

0.11 

-1.01 

0.1 

31 

0.73T0.21 

0.71T0.28 

0.98L0.55 

1.34L0.02 

39.7 

Ruprecht 79 

277.1 

0.08 

-0.82 

0.09 

53 

0.77T0.15 

0.72T0.27 

0.93L0.45 

1.51L0.04 

53.6 

Stock 13 

290.5 

0.06 

1.59 

0.07 

9 

0.71T0.17 

0.79T0.44 

1.12L0.63 

1.46L0.04 

38.9 

Stock 16 

306.15 

0.09 

0.06 

0.08 

19 

0.76T0.16 

0.73T0.32 

0.96L0.52 

1.41L0.05 

43.9 

Trumpler 18 

290.98 

0.09 

-0.14 

0.09 

86 

0.73T0.19 

0.67T0.19 

0.91L0.40 

1.24L0.10 

65.1 

Trumpler 28 

356.01 

0.07 

-0.31 

0.07 

17 

0.81T0.11 

0.87T0.55 

1.06L0.69 

1.31L0.04 

18.5 

Berkeley 86 

76.66 

0.08 

1.28 

0.08 

48 

0.79T0.14 

0.73T0.30 

0.92T0.47 

1.31L0.04 


NGC 2244 

206.31 

0.25 

-2.06 

0.24 

256 

0.74T0.14 

0.68T0.13 

0.91L0.27 

1.39L0.05 


NGC 2264 

202.92 

0.36 

2.18 

0.33 

225 

0.71T0.16 

0.66T0.14 

0.93L0.30 

1.49L0.03 


NGC 6530 

6.08 

0.13 

-1.32 

0.12 

55 

0.87T0.10 

0.66T0.20 

0.77L0.28 

1.35L0.03 



Notes: The center of the spatial distribution of the cluster members is indicated in Columns 2-5 showing average galactic coordinates 
and respective 2cr deviations (A l and A b), which roughly indicate the cluster radius. 


For comparison among different samples, in Fig. [3 the 
dynamical time ( Td yn ) giv en by PZ MG10 was a dopt ed, as¬ 
suming that Ter = 2.8Tdyn dGieles fc Portegies Zwartll201ll b 
The Milky Way stellar groups are presented in two distinct 
distributions: Associations (II < 1) and YMCs (II > 1). 
In general, the objects in other Local Group galaxies, like 
LMC, SMC and M31, as well as outside the Local Group, 
seem to follow the same distribution as YMCs of Milky Way, 
but are more scattered. 

Figure [7]clearly shows that all clusters of the sample are 
more likely associations, following the same trend of Milky 
Way objects having II < 1 that corresponds to unbound 
expanding groups. No differences are found when comparing 
II with the clustering characteristics of our sample, meaning 
that both fractal and smooth radial density profile types of 
objects follow the same age vs. crossing time correlation. 

Among the clusters studied by 

ICartwright fe Whitworthl (|2004l b Taurus (II = 0.1) co¬ 
incides with the distribution shown by our sample and 
other associations, while p Oph (II = 1.5) and Chamaeleon 
(II = 0.74) are closer to IC 348, which is in the transition 
region (II = 1). IC 2391 (II = 21.2) is found in the region 
occupied by older clusters of other galaxies. The scattering 
observed in this region is probably due to observational 
difficulties to determine accurate effective radius of distant 
objects. Our sample is not affected by this observational 
constraint. 


5.3 Clustering Evolution 

The values of Q and age, observed in the sample, 
were compared with diff erent YSOs classes analysed by 
ISchmeia fe Klesser] d2006lb We noted that Q values of our 
clusters are systematically above those obtained with 3D 
calculations (Q ~ 0.7) for the YSOs in the same range of 
ages (2 to 7 Myr). ISchmeia fe Klesseiil d2006l . see their Fig. 
5) suggest that 3D-Q values tend to be lower than the pro¬ 
jected 2D calculations. Indeed, a better agreement is found 
by comparing the results with t he calculations for the pro - 
jections in 2D plans presented bv ISchmeia fc Klessenl d2006lb 

Considering that the evolution of clustering structure 
may be affected by the initial star-forming conditions, in 
Fig . [8] the resu lt s are compared with the calculations from 
IParker fc Dalel d2013l . see their Fig. 3) that used models 
with- and without feedback to simulate the effects of the pres- 
ence (or not) o f ionizing so urces affect ing the star formation 
dPate, Ercolano fc Bonnehl 120121 . l2013t l . Even taking in ac¬ 
count the large error bars on Q, as illustrated in Fig. U 
the distribution of Q values in the sample, as a function of 
age, is consistent with feedback models, in particular those 
assu ming low initial densities (Ei = 0.6 to 83 stars pc -2 ). 
IParker fc Dalel (|2013h suggest that the lowest density models 
are the mostly affected by feedback, which are the oniy ones 
retaining some substructure for 5 to 10 Myr. Since clusters 
with low densities have longer relaxation times, their sub¬ 
structures are not erased due to negligible stellar mixing. 

Other young clusters from literatu re were a lso in clude d 
in Fig. [8] Excepting Cygnus OB2 dWright et al.l 120141 1 . 
the other clusters like p Oph, Taurus, Chamaeleon, IC348 
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Figure 7. Distribution of age vs. crossing time. Th e sam ple is 
compared to results from I Cartwright Sz Whitworthl (120041) (☆) 
and from PZMG10 (same symbols as Fig. E}. The groups of Milky 
Way show two distributions: bound YMCs (+) and unbound as¬ 
sociations (x), separated by a dotted line (II = 1). 


TT >' 


D n m ^ 


** * * 

* ** 


< 


LD 



Age [Myr] 

Figure 8. Distribution of Q as a function of age for the sam¬ 
ple and other clusters (same as Fig. Cygnus OB2 is indicated 
with error bars. The results are compared to the simulations with 
feedback (solid lines) and without feedback (dashed lines), indi¬ 
cating the initial densitie s (number of stars pe r pc 2 ) for some of 
the models adapted from lParker fe Dalel d2013h ■ 


llSchmeia fc Klessenll2006l ). and a Ori (ICaballerol 120081 ) are 
similar to o ur sam ple, coinc iding with low de nsity feedback 
models from lDale, Ercolano fc Bonnelll d2012i . 120131 ). 

T he results are also comparable to those from 
IParker et al.l (|2014h for the evolution of Q based on simula¬ 
tions of different values of virial ratio a v i r , the total kinetic 
energy over the potential energy. Our clusters have a distri¬ 
bution very similar to the simulation of aw = 1.5 (see their 
Fig. 3g,h,i), which represents a globally supervirial fractal 
that is “hot” and unbound, according the terminology in¬ 
troduced by them. The same is found for Cygnus OB2, a 
very substructured cluster (Q = 0.4), which is better repre¬ 
sented by c alculations for unbound regions under supervirial 
conditions dWright et al.ll2014h . 


6 CONCLUSIONS 

In the present work, the methodology used by FGH12 was 
adopted to perform a fractal analysis of the sample stud¬ 
ied by SG12, aiming to enlarge our previous study for the 
clusters NCG 6530, Berkeley 86, NGC 2244, and NGC 2264. 
About half of the sample (52%) shows substructures, indi¬ 
cated by the measured parameter Q < 0.8, while the other 
48% of the o bjects are cen t rally concentrated (Q > 0.8). 
According to IParker et al.l (|2014j), regions having low Q 
(Q < 0.8 or 1) must be dynamically young. 

The projected spatial distribution of very likely mem¬ 
bers, which have membership probability P > 50%, shows 
a m vs. s correlation that is confirmed by results from lit¬ 
erature for other young clusters. On the other hand, the 
sub-sample of candidate members (P < 50%) do not follow 
this trend. 

Centrally concentrated clusters of the sample have s 
increasing with the mean spatial density (n), which is an 
expected result since s shows more variation for radial clus¬ 
tering. This could also be the explanation for the lack of 
trend in the n vs. s for objects with fractal clustering. How¬ 
ever, a Ori shows a too high surface density, being out of the 
correlation with s found by us. This discrepancy is possibly 
due to uncertainties on the parameters adopted for <7 Ori 
( Nt = 340 or R = 3 pc). 

Visual extinction maps of clouds nearby the clusters 
were used in order to study the geometric structure of gas 
and dust distribution, by measuring their fractal dimension 
D 2 ■ A transform function was adopted to compare the re¬ 
sults, measured on projected 2D maps, with simulations 
based on 3D models. Intermediate values 2.1 < D 3 < 2.5 
were found, which correspond more likely to substructures 
than to smooth radial profiles. 

The distribution of fractal dimension D 3 (estimated for 
the observed clouds) as a function of Q parameter (mea¬ 
sured for the clusters) shows that half of the sample tends 
to follow a theoretically e xpecte d Q x D 3 relatio n, cal- 
culated for artificial da ta llCartwright fc Whitworthl 120041 ; 
ISanchez fc Alfaroll2009l ). These clusters have Q < 0.78 that 
indicates the presence of substructures, similar to those ob¬ 
served in the nearby clouds. On the other hand, 12 clusters 
(48%), among them NGC 6530, have Q > 0.8, which corre¬ 
sponds to radial distribution of stars that does not coincide 
with the structure of the clouds. 

It should be expected that these two groups of clusters 
could have had different virial conditions in the early for¬ 
mation, which are related to different thermal paths: cold 
collapse for objects centrally concen trated , or warm c ol- 
lapse for those showing substructures llDelgado et al.ll2013i ). 
Therefore, differences on dynamical evolution could also be 
expected, but they were not observed for these two subsam¬ 
ples. 

The analysis of dependence of age on parameters like 
core radius and crossing time show that our clusters have 
many similarities among each other. They are also similar to 
other young stellar groups of Milky Way, in particular the 
unbound associations, as indicated by the dynamical ages 
n < 1 found for the entire sample. 

By comparing the evolution of Q with theoretical stud¬ 
ies, the clusters show a distribution similar to simulations 
that include feedback of ionizing source in the models 
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adopted bv lParker fe Dale! J2013h . According these authors, 
the calculations using low values of density are the mostly 
affected by feedback, corresponding to the only simulations 
retaining some substructure for 5 to 10 Myr, which appears 
to be the case of our sample. Other young clusters from lit¬ 
erature also coincide with these models. A slight tendency 
is observed for objects that have smooth density profile 
( Q > 0.8), like IC 348 and p Oph, being more similar to 
the calculation with Si = 24 stars pc -2 . On the other side, 
Taurus, which shows fractal structure, coincides with the Si 
= 0.6 stars pc -2 model. Our objects are found in between 
this range of densities. 

The observed distribution of Q as a function of age was 
also compared to the simulations bv lParker et alJ d2014|j as¬ 
suming different values of virial ratio ( a v i r )■ The results are 
quite similar to the calculations that use aw = 1.5, corre¬ 
sponding to “supervirial” fractal. 

Since no difference concerning the Q evolution is noted 
for our objects, independently if they are fractal or have 
radial density profile, it is not evident to state if they had 
or not similar conditions in the early formation. 

Additional discussions on the cluster formation condi¬ 
tions must include the initial density distribution of parental 
clouds, as well as t he effects of turbulence. According 
iGirichidis et alJ (i2012l h compressive modes in a flat density 
profiles tend to form substructures, instead of centrally con¬ 
centrated distributions of stars. The fractal dimension cal¬ 
culated for clouds, projected in the direction of the studied 
clusters, allow us to infer the geometric distribution of the 
dense regions. A more detailed study on dynamical condi¬ 
tions of the gas is required to better understand the scenario 
of clusters formation and their relation with their parental 
clouds. 
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